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A variety of observations indicate that the universe is dominated by dark energy with negative 
pressure, one possibihty for which is a cosmological constant. If the dark energy is a cosmological 
constant, a fundamental question is: Why has it become relevant at so late an epoch, making today 
the only time in the history of the universe at which the cosmological constant is of order the 
ambient density. We explore an answer to this question drawing on ideas from unimodular gravity, 
which predicts fluctuations in the cosmological constant, and causal set theory, which predicts the 
magnitude of these fluctuations. The resulting ansatz yields a fluctuating cosmological "constant" 

■ which is always of order the ambient density. 
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_ ■ I. INTRODUCTION 

(D ■ 

^/J • The most startling discovery to emerge from the recent plethora of cosmological data is that the universe appears 
, to be dominated by dark energy ^ ^ . We know that this dark energy accounts for roughly seventy percent of the 
I ' energy density in the universe, does not cluster like ordinary matter, and has negative pressure. Otherwise, we are in 
' the dark about the nature of this extraordinary phenomenon. 
T-H Perhaps the most popular explanation is that the dark energy is due to a cosmological constant, for such a parameter 

^ was introduced into general relativity at its birth and has remained an important tool for cosmologists seeking 
, to model the observed universe ||]. The strongest argument against the cosmological constant is that naively we 
' expect it to contribute an energy density, p\, of order — (87rG)~^, where G is Newton's constant and rUp is the 
: reduced Planck Mass^ . This estimate is some one hundred and twenty orders of magnitude larger than the observed 
value. An equivalent way of stating the problem is to note that only today is the cosmological constant of order the 
ambient density in matter or radiation. At all past epochs, pA was sub-dominant and immensely so. Many people 

■ have felt that no theory could naturally predict such a tiny value for A (or equivalently such a late epoch for it to 
(-H ' become relevant) without predicting A to vanish entirely, and for this reason they have sought other explanations of 
Q_i^ the observations. 

Many alternatives to the cosmological constant have been proposed. Most significant among these arc quintessence 
models in which the dark energy is due to a homogeneous scalar field shifted away from the true minimum of its 
potential Like a simple cosmological constant, many of these suffer from the "Why Now?" problem: Why does 
^ the quintessence field come to dominate only recently? They also typically need to explain the small mass scale 
ILJ necessary for the field to be important today (m < 10"^'^ eV). Even more disturbing, none of them are connected to 
. ,—{ [ realistic particle physics models. Perhaps then, instead of altering the energy content of the universe, we need to look 
' in another direction and modify gravity in order to explain the dark energy today. 
5h I The simulations reported here flesh out an old heuristic prediction of a fluctuating cosmological term arising 
from the basic tenets of causal set theory.^ In normal usage, the words "cosmological term" refer to a contribution 
to the effective stress-energy-momentum tensor of the form T^^ = (7^yA(x). However, in classical General Relativity 
(GR) such a A(a;) must be constant if the total energy momentum in other matter components is separately conserved. 
Here we consider a specific modification of GR motivated by the search for a theory of quantum gravity based on 
causal sets. 

Although the ultimate status and precise interpretation of the prediction of a fluctuating A must await the devel- 
opment of a quantum dynamics for causal sets , the basic lines of the argument are simple and general enough that 
they have a certain independence of their own. In this paper we review the motivation for a fluctuating A from causal 
set theory, propose an ansatz for the form of these fluctuations, apply the latter to the Friedmann equation with 
time-dependent cosmological term, and find that we can have a viable cosmology for some fraction of the solutions. 
Finally, we address issues related to our choice of evolution equations. 
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^ In this paper, we will use units in which h = c = rap 
^ For an introduction to the causal set hypothesis see 
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II. CAUSET THEORY 

Here, we will review the arguments leading to a fluctuating cosmological term and then describe the specific ansatz 
via which we have chosen to implement their main implication: that A can be expected to fluctuate and with a 
magnitude that diminishes as the universe grows older. 

In causal set ("causet") theory, the predicted fluctuations arise, as a kind of residual (and nonlocal) quantum 
effect, from the underlying space-time discreteness. More specifically, the basic inputs to the argument are: space- 
time discreteness leading to a finite number N of elements; the interpretation of space-time volume V as a direct 
reflection of N; the conjugacy of A to space-time volume V; and the existence of fluctuations in V coming from 
Poisson fluctuations in N. (Of these four inputs, the first is not peculiar to causal sets, but the remaining ones all 
are to a greater or lesser extent.) 

The two most basic tenets of causal set theory are first, that the causal ordering of events in macroscopic space-time 
reflects a more fundamental order relation among the elements of an underlying discrete structure to which continuous 
spacetime is only an approximation, and second, that the four-volume of a region of spacetime reflects the number of 
discrete elements of which the region is "composed". The hypothesized discrete substratum or causal set is taken to be 
a partially ordered set and its dynamics is conceived of as a kind of growth process in which elements come into being, 
one at a time. Although a classically stochastic dynamics expressing these ideas is by now fairly well developed , 
a corresponding quantum dynamics is only just beginning to be sought. Any prediction of quantum fluctuations in A 
must therefore rest on an anticipation of certain features of this "new QCD" (quantum causet dynamics). 

Let us begin by assuming that, at some level of approximation, this dynamics will correspond to a space-time "path 
integral" in which one is summing over certain classes of four-geometries. At the deeper level however, this will still be 
a sum over causal sets. Then let us take from the already developed classical growth models for causal sets the feature 
that the ever growing number N of causet elements plays the role of a kind of parameter time - the time in which the 
stochastic process which mathematically represents the growth unfolds, and with respect to which the probabilities 
are normalized. Just as one does not sum over time in ordinary quantum mechanics, one would not expect to sum 
over causets with different values of N in the quantum theory. But, because number corresponds macroscopically to 
volume V, this translates into the statement that one should hold V flxed in performing the gravitational path integral. 
Any wave function that arises will therefore depend not only on suitable boundary data (say a three-geometry) but 
also on a four- volume parameter V. Such a restricted path integral may be called "unimodular" . 

Now the unimodular modiflcation of ordinary GR has been fairly well studied jlll , and it is understood that within 
it, A and V are conjugate in the same way that energy and time are conjugate in ordinary quantum mechanics. (Indeed 
this is almost obvious from the fact that the cosmological constant term in the action-integral of general relativity is 
just the product — AV.) In particular, this means that, to the extent that V is held fixed in the gravitational sum- 
over-histories, A will be entirely undetermined by the fundamental parameters of the theory. (Again this is almost 
obvious by reference to the classical limit of unimodular gravity, where the Lagrange multiplier used to implement 
the flxed V constraint combines with any "bare" A in such a way that the observed or "renormalized" A represents 
nothing more than a constant of integration.) 

If this were the whole story, then our conclusion would be that A is subject to quantum fluctuations (just like 
energy E in ordinary quantum mechanics) but it would not be possible to say anything about their magnitude, nor 
about the magnitude of the mean A about which the fluctuations would occur. 

But here there enters a second aspect of the causal set hypothesis that we have not mentioned earlier. In order to 
do justice to local Lorentz invariance, the correspondence between number and volume cannot be exact, but it must 
be subject to Poisson type fluctuations^, which of course have a typical scale of ^/N. This means that, in holding 
N fixed at the fundamental level, we in effect fix V only up to fluctuations of magnitude ±VV. (Notice that these 
are not dynamical fluctuations. Rather they occur at a kinematic level: that of the correspondence between order 
theoretic and spatio-temporal variables.) Hence, we do end up integrating over some limited range of V after all, and 
correspondingly we do determine A to some degree — but only modulo fluctuations that get smaller as V gets larger. 
Speciflcally, we have 

AA- 1/AV- l/\/v. (1) 

As any proper dilemma should, that of the cosmological constant has two horns: Why is A so nearly zero and Why 
is it not exactly zero? None of what we have said so far bears on the first question, only on the second. All we can 



More specifically, the correspondence between the underlying causet and the approximating space-time is via a notion of "Poisson 
sprinkling" at unit density, see references M, IVM for details. 
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conclude is that partially integrating over V in the effective gravitational path integral will drive us toward some value 
of A. We must assume, as the evidence overwhelmingly suggests, that this "target value" is zero, for reasons still to 
be understood.'' Then we end up predicting fluctuations about zero of a magnitude given by 

Independent of specifics, the space-time volume V should be roughly equal to the fourth power of the Hubble radius, 
H~^. Therefore, at all times we expect the energy density in the cosmological constant to be of order 

the critical density (recall that we are setting 8ttG = 1). We thus obtain a prediction for today's A which agrees in 
order of magnitude with current fits to the astronomical data. And this argument is not limited to today: at all times 
we expect the energy density in the cosmological constant to be of order the critical density. 

This is the basic idea, but any attempt to implement it immediately raises questions whose answers we can at 
present only guess at, pending the development of a fuller quantum dynamics for causets. At a conceptual level there 
is first of all the question of precisely how to interpret the V that figures in equation (|^) and second of all the question 
how to incorporate a fluctuating A into some suitable modification of the Einstein equations. At a more practical 
level, if we aim to understand, for example, how fluctuations in A would have affected structure formation, we need to 
know, not only their typical magnitude at each moment of cosmic time, but also how the fluctuations at one moment 
correlate with those at other moments. 

Concerning the conceptual questions, we will, for present purposes, resolve them provisionally as follows. First we 
impose spatial homogeneity, so that the Einstein equations reduce to a pair of ordinary differential equations for the 
scale factor a. Of these two equations, one, the so called Friedmann equation or Hamiltonian constraint^ is first order 
in time and embodies the energy law in this setting, while the second involves d and, in the case of a non-fluctuating 
A, adds no information to the first, except at moments when d = 0. They cannot both be compatible with a time 
dependent cosmological term when other energy momentum components are separately conserved, so we choose one 
over the other. Specifically, we choose to interpret A via its role in the Friedmann equation. That is, we retain the 
Friedmann equation but let A be a function of time, dropping the second equation entirely. 

The quantity V which governs the magnitude of the fluctuations in A we will identify (up to an unknown factor 
of order unity) with the volume of the past light cone of any representative point on the hypersurface for which we 
want the value of V, as illustrated in Figure |l|. Although this interpretation is somewhat at odds with the meaning 
that V has in the unimodular context, it seems more in accord with causality, and it is the only number accessible to 
observation in any useful sense. 

With these choices made, the only remaining question is what sort of random process we want to use to simulate our 
fluctuating A. Ideally perhaps, this would be some sort of "quantal stochastic process" (since the underlying process 
is quantal), but here we do the simplest thing possible and let the fluctuations in A be driven by those of an unadorned 
random walk. In fact the ansatz we will use has some appeal in its own right as an independent "story" of why the 
cosmological constant might be expected to fluctuate in any discrete quantum gravity theory that incorporates the 
equality N = V between volume and number of elements. 

With reference to the Einstcin-Hilbert Lagrangian, one could describe the cosmological constant as the "action 
per unit space-time volume which is due just to the existence of space-time as such, independent of any excitations 
such as matter or gravitational waves" . Re-interpreting volume as number of elements, we can say then that A is 
the "action per element". One would expect this to be of order unity in fundamental units, and if we identify the 
latter with Planck units, we get the old answer which is off by some 120 orders of magnitude. On the other hand, if 
we suppose that each element makes its own contribution and these contributions fluctuate in sign^ then the relative 
smallness of A will be explained; but one would also expect a residual ^/N contribution to S to remain uncanceled. 
Consequently, there would remain a residual contribution to the action per element of \/N /N = 1/ ^/N, in agreement 
with our earlier argument. 

To implement such an ansatz is now straightforward. What we need for the sake of the Friedmann equation is just 
A as a function of N (or equivalently of V). To produce such a function we just generate a string of random numbers 
of mean and standard deviation 1 (say) and identify A(A^) with the ratio S{N) /N, where S{N) is the sum of the 
first N of our random numbers. Modulo implementational details this is the scheme we have used in the simulations 
on which we report next. 



* One possible mechanism is that only A = is stable against the destructive interference induced by non-manifold fluctuations of the 
causal set. 

^ It would probably be more suitable to speak not in terms of action S but rather exp(iS'/?i) and say that the contributions (now 
multiplicative rather than additive) fluctuate in phase. 
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FIG. 1: Schematic representation of the backward light-cone at two difTerent cosmic times. Evolution of the scale factor between 
the two time slices is determined by the Friedmann equation while A varies stochastically. 



III. SIMULATIONS 



We take as the space-time volume 



m = 



4n 



i\3 



dt'a[t') 



t -|3 

dt"/a{t") 



(3) 



where a{t) is the scale factor of the universe at proper time t. Note from this formula that the backward light-cones 
depicted in Figure |l| are quite deceptive: because a(t) was much smaller in the past and vanishes at the Big Bang, 
most of the four- volume V of these light cones accumulates recently. One consequence of this is that V ^ H^'^ recently 
even if there was a period of cosmic inflation in the early universe. 

Our algorithm for calculating the cosmological constant at time-step i -f 1 is then to set 



6N, = N,+i - TV, - V(i,+i) - V{h) 



(4) 
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and then write 




(5) 



Here a is an unknown dimcnsionless parameter which governs the dynamics of the theory; ^^^i is a random number 
with mean and standard deviation 1; and 5*0 is set to zero at some very early time Iq. We then expand the universe 
according to 



recompute the new space-time volume and repeat. 

Figure |^ shows the evolution of the energy density in one such realization. During the radiation era, pA scales 
roughly as a""*, while during the matter era it scales as . Thus at all times it is comparable to the ambient 
energy density. If the recipe we have devised for implementing the ideas of causal set theory and unimodular gravity 
is an accurate approximation to the ultimate quantum theory, then these modifications of GR do indeed lead to an 
Everpresent A, a cosmological term which is always with us [[l3|. 



FIG. 2: Evolution of the energy densities in the universe. The thick curve is the absolute value of the energy density in the 
cosmology constant. The fluctuating pA is always of order the ambient density, be it radiation (early on) or matter (later). 
Here the dimensionless parameter a which governs the amplitude of the fluctuations has been set to 0.01. 

Hidden in the gross structure of Figure ^ are the fluctuations about this average scaling. These fluctuations are 
crucial if the theory is to describe the real universe for two reasons: First, there cannot be too much excess energy at 
a ^ 10~^ or else the successful predictions of Big Bang Nucleosynthesis (BBN) will be destroyed. Second, if pA scales 
exactly as matter today, it will not have the correct equation of state to account for the cosmological observations. 
Figure ^ shows the ratio of the energy density in A to the total energy density as a function of the scale factor for 
another realization, this time with a slightly larger value of a. This ratio, Q\, fluctuates about zero with an amplitude 
of order 0.5 (as we will shortly see, this amplitude is a function of a). In this particular realization, A accounts for 
over fifty percent of the energy density today and changes very little going back to redshift z — 1 (a = 0.5); thus it 
behaves recently as a true cosmological constant, and therefore satisfies the observed cosmological constraints. 

In half the realizations, p\ will be positive today. Whether or not it is positive enough to explain the observations 
then becomes a question of probability. For a = 0.02, it clearly is not that improbable (indeed in the same run, we 
see another spike in the energy density at a ~ 0.1). 




(6) 
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FIG. 3: The ratio of the energy density in cosmological constant to the total density as a function of scale factor. Here a — 0.02. 

The same qualitative argument applies to the BBN constraint. In fact the situation there is even better. Half the 
time the extra energy density will be negative, thereby reducing the total energy density in the universe. This in 
turn will reduce the predicted abundance of *He. There is some disagreement at present as to whether the current 
observations agree with the standard cosmological model or not , with some cosmologists arguing that the observed 
abundances are too low. A negative pA fixes this problem. 

Why have we chosen a to be small? Our choice is in response to a fundamental incompleteness in our implemen- 
tation. If a « 1, there will inevitably be times during which the total effective energy density, the sum of the terms 
in parenthesis on the right side of equation (|^), goes negative, thereby invalidating the equation. (Whenever this 
happens, we terminate the run.) In the next section we offer some thoughts on this problem; here we simply spell it 
out. 

Figure || shows a history for a = 0.02 going back to the time of decoupling. If we had started earlier, say at the 
"Planck time" a = 10~^^, we would have had only about a 1 in 3 chance of completing the run without hitting a time 
at which ptot went negative. Moving a down to 0.01 evades this problem; for that parameter choice, very few runs hit 
a time at which the total energy goes negative. However, for a that small, the fluctuations are also small. Figure ^ 
shows a histogram of final values of ^Ia for 6000 realizations each with a = 0.01. Only rarely does the final value of 
f^A approach those necessary to explain the observations. 

There is therefore a tension: if we push a too low, it becomes very unlikely that pA will be large enough today to 
agree with observations. If we push a too high, there inevitably comes a time at which the total energy density in the 
universe becomes negative, and the simulation cannot continue. Of course we are dealing with probabilities, so for 
any value of a there is always the chance that the total energy density remains positive throughout the history of the 
universe and the final value of pA is large enough to account for observations. Fortunately, this happens reasonably 
often for a in the range 0.01 — 0.02. Nonetheless, we suspect that we will ultimately have to deal more directly with 
possibility that ptot goes negative. 

IV. COMPLICATIONS 

We can think of two ways to deal with a negative ptot without having to terminate the simulation: change the 
implementation so that this never occurs or reinterpret ptot going to zero (or negative) so as to give a viable cosmology 
without having to fine tune a. 

One approach would be to suppose that A fluctuates but is positive semi-definite. This is the position adopted by 
Ng and van Dam in reference [[l5|. There, they argue that the kernel for the Euclidean gravitational path integral 
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FIG. 4: A histogram of the final value of D,a, the ratio of pA to the total density. The dimensionless parameter governing the 
fiuctuations in A has been set to a = 0.01. 



over A histories takes the form 

e"'^^ cx exp ( — ^ 1 , (7) 

in our units. From this, they argue that the most probable value of A is zero and that, if it is not zero, it must be 
positive. As they observe, however, this result is peculiar to the assumptions of Euclidean quantum gravity with all its 
uncertainties and controversy. In particular, this result does not seem to follow from causal set theory or unimodular 
gravity by themselves and we do not favour it. 

Another possibility would be to suppose that the cosmological term comes from a decrease in the local energy 
of one of the matter fields or gravitational waves. This is the philosophy of, for example, Chen and Wu who 
consider a non-fluctuating, but time dependent, cosmological term A(t) cx a^^(i). Of course, this supposition is forced 
upon us if all of Einstein's equations are to be simultaneously satisfied exactly. That is, Einstein's equations — the 
contracted Bianchi identity in particular — require that total energy-momentum be conserved. Thus, in classical GR, 
the cosmological term cannot fluctuate without a compensating fluctuation in the energy-momentum density of one 
or more of the matter fields. 

As stated in the introduction, our approach has been to solve this problem by maintaining only the Friedmann 
equation as exact. Nevertheless, let us try instead to adapt the solution above to our case. Suppose that some matter 
component, let us take gravitational waves as a concrete example, is somehow converted into the energy density of a 
cosmological term, while the energy density in every other component (dust, radiation, etc.) is separately covariantly 
conserved. The first law, ^ [(ptot +Ptot)a'^] = a^^^^j^, apphed to these two components becomes: 

Pgrav = -4i7pg,.av ^ A(t) , (8) 

where H is the Hubble parameter. As could be expected, an increase in A must lead to a decrease in the energy 
density in the gravity waves. However, for a generic fluctuating A(t), the cosmological term might increase enough 
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that the energy density in gravity waves becomes negative. It's not clear how this could be interpreted and it appears 
that we would simply exchange one problem with another. Note that, in the case of Chen and Wu this is not 
a problem. Their cosmological term decreases monotonically with the expansion of the universe. Thus, we see that 
this solution can work with a fluctuating A if we demand that A < 0. In fact, relaxation processes of this sort have 
been considered for some time. The earliest we are aware of is that of Abbott Recently, there has been renewed 
interest in a similar suggestion of Brown and Teitelboim where domains of four-form flux decay spontaneously, 
relaxing the effective local value of the cosmological term, see for recent references. The difficulty with these 
proposals is again the "Why now?" problem: relaxation rates and/or boundary values must be tuned for any hope 
to obtain a viable cosmology.^ 

All in all, neither of these proposals really seems to address the central difficulty: Within the contexts of causal set 
theory and unimodular gravity, the sign of the total effective energy density is fundamentally not constrained. We see 
no good reason to assume either that A is positive definite or that A will always decrease. Thus, let us seek instead to 
understand what happens when ptot approaches and perhaps goes through zero. Our guide will be the classical theory. 
Consider, for example, a dust filled, flat universe with a negative cosmological constant Ao < (a true constant). The 
0-0 component of Einstein's equations gives us the Friedmann equation: — ptot- Meanwhile, the i-i component 
gives us the deceleration: 2| = — |Ao| — . We see that, once the matter has red-shifted enough that the total energy 
density vanishes, the universe stops expanding and begins to contract. As it contracts, the energy density in matter 
once again begins to exceed the magnitude of the cosmological term and ptot never becomes negative. 

We expect that something like this phenomenology will carry over into our case except that, with a fluctuating 
cosmological term, it seems likely that the collapse can reverse itself if the cosmological term later becomes sufficiently 
negative a second time. This is in contrast with the classical example above where once the universe starts to collapse, 
the matter term always dominates and keeps a /a negative. In our model, however, we would expect the amplitude of 
the cosmological term's fluctuations to track the matter or radiation energy density in the collapsing universe. 

Of course, none of this follows from the simple evolution ansatz we have applied in this paper. Such detailed 
dynamical understanding must await further developments. Nevertheless, it is reasonable to suppose that the complete 
theory will reduce in stages: a full theory with non-metric structures at the Planck scale and a semi-classical theory 
describing metric structures at larger scales. Furthermore, it is reasonable that the semi-classical theory will be 
describable as some sort of sum-over-histories, where the intermediate states are three-geometries. We can envisage 
the evolution equation we propose as some sort of classical approximation to propagation in a stochastic potential. In 
the sum-over-histories theory, we expect the case ptot < will correspond to a tunneling-type solution. Our difficulty 
in handling ptot < here is, in this sense, no different from the standard problem of finding an effective, classical 
description of barrier penetration. 

V. CONCLUSION 

It is still too early to understand the full implications of recent cosmic discoveries that point to dark energy in 
the universe. A number of possibilities have previously been explored in detail, including a non-zero cosmological 
constant A and zero A with dark energy hidden in a scalar field. 

It is also possible, though, that the measurements are telling us that we need to modify our understanding of space 
and time. In particular, the notion that space-time is continuous may be simply an approximation that breaks down 
on scales as small as the Planck scale. If so, drawing on ideas from causal set theory - which postulates a discrete 
space-time - and unimodular gravity, we have shown that the cosmological "constant" need not be a fixed parameter. 
Rather, it arguably fluctuates about zero with a magnitude l/\/V, V being some measure of the past four volume. The 
amplitude of these fluctuations is then of the right order of magnitude to explain the dark energy in the universe. This 
argument is so general that it would apply at all times, and, indeed, we expect the energy density in the cosmological 
"constant" to always be of order the ambient density in the universe. 

In §IV, we presented a number of issues which inevitably will confront anyone wishing to implement this idea. Until 
these issues are resolved, it will be difficult to make unambiguous, robust predictions. Nevertheless, one can already 
see that this theory of a fluctuating A differs signiflcantly from most other solutions to the dark energy problem. 
Most important for its testability is the notion that it may have affected the evolution of the universe at early times. 
Thus, the primordial generation of perturbations during a possible inflationary phase; production of light elements in 



^ Of course, back when Abbott and Brown and Teitelboim first made their suggestions, there was no compelling evidence for a non- 
vanishing cosmological term. One needed only to make it small enough. 
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Big Bang Nucleosynthesis; acoustic oscillations in the background radiation; and the evolution of structure at more 
recent times all may yield clues and tests of the idea of an everpresent A. 
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